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Abstract. The partially asymmetric exclusion process (PASEP) is an impor- 
tant model from statistical mechanics which describes a system of interacting 
particles hopping left and right on a one-dimensional lattice of n sites. It is 
partially asymmetric in the sense that the probability of hopping left is q times 
the probability of hopping right. Additionally, particles may enter from the 
left with probability a and exit from the right with probability /3. 

In this paper we prove a close connection between the PASEP and the 
combinatorics of permutation tableaux. (These tableaux come indirectly from 
the totally nonnegative part of the Grassmannian, via work of Postnikov, and 
were studied in a paper of Steingrimsson and the second author.) Namely, 
we prove that in the long time limit, the probability that the PASEP is in a 
particular configuration r is essentially the generating function for permutation 
tableaux of shape A(r) enumerated according to three statistics. The proof of 
this result uses a result of Derrida, Evans, Hakim, and Pasquier on the matrix 
ansatz for the PASEP model. 

As an application, we prove some monotonicity results for the PASEP. We 
also derive some enumerative consequences for permutations enumerated ac- 
cording to various statistics such as weak excedence set, descent set, crossings, 
and occurences of generalized patterns. 



1. Introduction 

The partially asymmetric exclusion process (PASEP) is an important model 
from statistical mechanics which is quite simple but surprisingly rich: it exhibits 
boundary-induced phase transitions, spontaneous symmetry breaking, and phase 
separation. The PASEP is regarded as a primitive model for biopolymerization 
[llj . traffic flow [T5], and formation of shocks [5]; it also appears in a kind of 
sequence alignment problem in computation biology [4]. More recently it has been 
noticed that the PASEP model has relations to orthogonal polynomials [13] , and to 
interesting combinatorial phenomena [§1 [5J [Tl] . The goal of this paper is to prove a 
precise connection bctwcccn the PASEP model and permutation tableaux, certain 
0—1 tableaux introduced in [TB] . 

The PASEP model describes a system of particles hopping left and right on a 
one-dimensional lattice of n sites. Particles may enter the system from the left with 



2000 Mathematics Subject Classification. Primary 05E10; Secondary 82B23, 60C05. 
Key words and phrases, permutation tableax, asymmetric exclusion process, matrix ansatz, 
Eulerian numbers, generalized patterns. 



2 



SYLVIE CORTEEL AND LAUREN K. WILLIAMS 



a rate adt and may exit the system from the right at a rate [3dt. The probability 
of hopping left is q times the probability of hopping right. 

Let denote the probability that the PASEP model is in a particular configu- 
ration t in the steady state. Here, Z n is the partition function of the PASEP model. 
In this paper we consider a solution {D\, Ei, Vi, W\) to the "matrix ansatz" for the 
PASEP model which naturally relates to permutation tableaux. That is, we find 
matrices D\ and E\ and vectors V\ and W\ which satisfy the relations of Theorem 
12.31 we then prove that expressions of the form Wi(Y\^ =1 (T i D 1 + (1 — Ti)Ei))Vi, 
where r = (ti, . . . ,t„) £ {0, 1}", are (Laurent) polynomials in q,a,(3 which enu- 
merate permutation tableaux of a fixed shape A(r) according to three statistics. 
Using a result of Derrida et al [7j, we are then able to prove our main result: that 
f T (q) is the generating function for permutation tableaux of a fixed shape A(r). 

It then follows from work of the second author and Steingrimsson [TO] that f T (q) 
is also the generating function for: permutations in S n +i with a fixed set W(t) of 
weak excedences, enumerated according to crossings; permutations in S n +i with a 
fixed set D(t) of descents, enumerated according to occurrences of the generalized 
pattern 2 — 31. Additionally, these results imply that the expression W\{D\+ E\) n Vi 
for the partition function Z n is also the weight-generating function for permutations 
in SVi+ij enumerated according to crossings, and permutations in S n +\, enumerated 
according to occurrences of the generalized pattern 2 — 31. 

Our main result refines the theorem of the first author [5], who showed that 
if a = (3 = 1, then in the steady state, the probability that the model is in a 
configuration with k occupied sites is equal to Ek + 1 -£+ 1 ^ ) w here Ek+i, n +i(<l) is the 
q-Eulerian polynomial introduced by the second author [17] . 

The structure of this paper is as follows. In Section [2] we define the PASEP 
model and review the "matrix ansatz," presenting a solution (D\, E\) which has an 
interpretation in terms of permutation tableaux. In Section [3] we define permutation 
tableaux, certain 0—1 tableaux which were defined in [16] and which are naturally in 
bijection with permutations. We then show how the solution (D\, E\) to the matrix 
ansatz leads to a natural connection between the PASEP model and permutation 
tableaux, and hence to permutations. In Section 01 we prove the main result of the 
previous section, and in Section [5l we give some applications. Finally, in section [6j 
we show how a classical solution (Dq,Eq) of the matrix ansatz leads to a natural 
connection between the PASEP model and bicolored Motzkin paths. This recovers 
results of Brak et al [3] ■ 

It is interesting to note that permutation tableaux are closely connected to total 
positivity for the Grassmannian |17[ 1 1 2 j . This suggests an intriguing connection 
between total positivity and the PASEP model. 

Acknowledgments: We would like to thank the referee for insightful com- 
ments. The second author is grateful to Persi Diaconis for interesting discussions 
in Banff and for encouraging her to work on this project. She would also like to 
thank Kolya Reshetikhin and Richard Stanley for encouragement and useful com- 
ments. 
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2. The PASEP model and the "matrix ansatz 



In the physics literature, the PASEP is defined as follows. 

Definition 2.1. We are given a one- dimensional lattice of N sites, such that each 
site i (1 < i < N) is either occupied by a particle fa = 1) or is empty fa = 0). At 
most one particle may occupy a given site. During each infinitesimal time interval 
dt, each particle in the system has a probability dt of jumping to the next site on its 
right (for particles on sites 1 < i < N — 1) and a probability qdt of jumping to the 
next site on its left (for particles on sites 2 < i < N ). Furthermore, a particle is 
added at site i = 1 with probability adt if site 1 is empty and a particle is removed 
from site N with probability [3dt if this site is occupied. 

Note that we will sometimes denote a state of the PASEP as a word in {0, 1} 
and sometimes as a word in {o, »} N . In the latter notation, the symbol o denotes 
the absence of a particle, which one can also think of as a white particle. 

It is not too hard to see [9] that our previous formulation of the PASEP is 
equivalent to the following discrete-time Markov chain. 

Definition 2.2. Let a, (3, and q be constants such that < a < 1, < /3 < 1, and 
< q < 1. Let Bn be the set of all 2 N words in the language {o, •}*. The PASEP 
is the Markov chain on E>n with transition probabilities: 

• If X — A • oB and Y — A o uB then Pxy = N X +1 (particle hops right) and 
Py.x = ]v+T (particle hops left). 

• If X = oB and Y = uB then Px.y = ivTT (P ar ticle enters from left). 

• If X — B* and Y = Bo then Px.y = N ^ +1 (particle exits to the right). 

• Otherwise Px,y = for Y ^ X and Px.x = 1 — J2x^y Px,y- 

See Figure [T] for an illustration of the four states, with transition probabilities, 
for the case n = 2. Note that the probabilities of the loops are determined from 
the figure by the fact that the sum of the probabilities on all outgoing arrows from 
a given state must be 1. So for example the probability on the bottom- most loop 

j s 1 _2 _ a _ §_ _ 3-q-a-P 
IB x 3 3 3 3 



In the long time limit, the system reaches a steady state where all the proba- 
bilities P n (ri, T2, . . . ,T n ) of finding the system in configurations (ti,T2, . . . , r n ) are 




o 
•o 



coo 



• • 



om 



Figure 1. The state diagram of the PASEP model for n — 2 
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stationary, i.e. satisfy 

d 

— P„(ri, ...,T n ) = 0. 

Moreover, the stationary distribution is unique [TJ, as shown by Derrida et al. 

The question is now to solve for the probabilities P n (n, . . . , r„). For convenience, 
we define unnormalized weights f n (Ti, ■ ■ ■ , t„), which are equal to the P n {j\, . . . , T n ) 
up to a constant: 

P„(ti, . . . ,t„) = / n (n, . . . ,T n )/Z n , 

where Z n is the partition function ^ r f n {Ti, ■ ■ . , t„). The sum defining Z n is over 
all possible configurations r £ {0, 1}™. 

The "matrix ansatz" has been used by Derrida et al [7J to obtain exact expres- 
sions for all the P n (ti , . . . , r n ) . 

More precisely, they show the following. 

Theorem 2.3. [7 j Suppose that D and E are matrices, V is a column vector, and 
W is a row vector, such that the following conditions hold: 

DE - qED = D + E 
DV=hf 

WE = -W 

a 

Then 

n 

fn(n, .. . ,t„) = WiHinD + (l - n)E))V. 

i=l 

Note that Y[i=i( T i^ + (1 — T i)^) i s simply a product of n matrices D or E with 
matrix D at position i if site i is empty (r^ = 0). 

Remark 2.4. It follows from Theorem \2.3\ that the partition function Z n is equal 
to W{D + E) n V. 

We will now describe a solution {Di,E\, V\, Wi) to the matrix ansatz. It seems 
to have not been considered before, although it naturally generalizes a solution 
given in [7] for the q = case. Our solution has an interpretation in terms of 
permutation tableaux and we found it through consideration of the combinatorics 
of these tableaux; however, one could also arrive at these matrices by choosing the 
basis {Vi, EVi, E 2 Vi, . . . } for the infinite-dimensional vector space that D and E 
act on, and then writing D and E according to this basis. (This method actually 
yields the transposes of the matrices D\ and E\ which we use below.) 

Let D\ be the (infinite) upper triangular matrix (d^) such that d^i+i = 
and d i: j = for j ^ i + 1. 
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That is, D\ is the matrix 
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Let Ei be the (infinite) lower triangular matrix (e^) such that for j < i, = 

^ (a" V" 1 (pi) + Eill ti +r h r )- Otherwise, e l3 = 0. 
That is, Ei is the matrix 



/ a- 1 

a- 1 /3 l + a- 1 q 

a" 1 /? 2 /3(l + 2a- 1 g) 1 + <? + a" V 

a" 1 /? 3 ^(l + Sa- 1 ^) + 2g + 3aT V) l + g + g 2 + arV 

V ; : ; ; 



/ 



Observe that when a = (3 = 1, we have 



7i_j)i ■ Here, [i]^ represents 
the fcth derivative of [i] with respect to q, and [i] is the g-analog of the number i, 
namely 1 + q + ■ ■ ■ + q z ~ 1 . 



And then Ei becomes the matrix 
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Let Wi be the (row) vector (1, 0, 0, ... ) and V\ be the (column) vector (1, 1,1,...). 
It is now easy to check that the required relations hold. 

Lemma 2.5. With the definitions of Di, Ei,Vi,Wi above, the following relations 
hold: DiEi - qEiDi = D 1 + E 1; DVi = \V X , and W X E = \W X . 



Proof. First note that 



W- 1 



i-2 

£ 

r=0 



when j < i + 1, and is equal to otherwise. Then note that 
q^DJij = g/T 1 (£!);,,_! = q ^"V" 2 
when 1 < j — 1 <i, and is equal to otherwise. 



"- + ^ 

J 7 r=0 



i — j + r + 1 
r 



i — j + r + 1 
r 
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Putting these together, we find that (D\Ei — qEiD\)i,j is equal to 



q 




v- 1 : + ^ E r r ^ s v 

V / r=0 V ' / s=l V ' 



7 — i / z — ' \ r 



when 1 < j < i. Also, we have that (D\Ei — qE\D\)i t j = /3 _1 when j = i + 1, is 
equal to a -1 /? 1-1 when j = 1, and is equal to when j > i + 1. These are precisely 
the matrix entries of D\ + E\, □ 

As we will show in Sections [3] and 2] the matrix product W / i(J|™ =1 (TiZ3i + (1 — 
Ti)Ei))V\ from Theorem 12.31 has a combinatorial interpretation as a generating 
function for permutation tableaux. 

Example 2.6. Suppose that we are considering the PASEP model with 3 sites, 
and are interested in computing the probability that in the long time limit, the sys- 
tem reaches the configuration which has a black particle in the second site (but the 
first and third sites are empty). This configuration is represented by the vector 
t = (0,1,0). Using Theorem \2.3[ we see that /3(0,1,0) is equal to the matrix 
product W\{E\DxE{)Vi; this product is cT 2 + a -1 /? -1 + a~ 2 (3~ 1 q. And by Re- 
mark [^m the partition function Z% is equal to W\(D\ + E%) V\, which is in this 
case a- 3 + 2cT 2 + 2a' 1 + a" 2 ^ 1 + 2a" 1 + 2/3" 1 + or 1 ft- 2 + 2/T 2 + f3~ 3 + 
q{cT 2 + oT 2 ^- 1 + 4ar 1 /T 1 + cT 1 ^ 2 + ft- 2 ) + q 2 ( a - 2 l3~ l + a^P' 2 ). Therefore 
the probability that in the long time limit the system is in the configuration (0, 1, 0) 
is 



■- 2 + a - 1 l3- 1 + a - 2 p- 1 q 



3. Connection with permutation tableaux 

Recall that a partition A = (Ai, . . . , A&) is a weakly decreasing sequence of non- 
negative integers. For a partition A, where A; = m, the Young diagram Y\ of 
shape A is a left-justified diagram of m boxes, with Ai boxes in the i-th row. We 
define the expanse of A or Y\ to be the sum of the number of rows and the number 
of columns. Note that we will allow a row to have length 0- i.e. we allow A^ = 
and we distinguish two partitions that differ in their number of empty rows. 

We will often identify a Young diagram Y\ with expanse n with the lattice path 
p(A) of length n which takes unit steps south and west, beginning at the north-east 
corner of Y\ and ending at the south-west corner. Note that such a lattice path 
always begins with a step south. See Figure [2] for the path corresponding to the 
Young diagram of shape (2, 1, 0). 
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If t € {0, l}"" 1 , then we associate to it a Young diagram A(t) with expanse n as 
follows. First we define a path p = {pi,.. . ,p n ) £ {S, W} n such that p\ = S, and 
Pi+i = S if and only if Tj = 1. We then define A(r) to be the partition associated 
to this path p. We denote by 4> the inverse of the above map: it is a bijection from 
the set of Young diagrams with expanse n to the set of n — 1-tuples in {0, 

As in [16) , we define a permutation tableau T to be a partition A together with 
a filling of the boxes of Y\ with 0's and l's such that the following properties hold: 

(1) Each column of the rectangle contains at least one 1. 

(2) There is no which has a 1 above it in the same column and a 1 to its left 
in the same row. 

We call such a filling a valid filling of Y\ . 

Note that the second requirement above can be rephrased in the following way. 
Read the columns of a permutation tableau T from right to left. If in any column 
we have a which lies beneath some 1, then all entries to the left of (which are 
in the same row) must also be 0's. 

Note that if we forget the requirement (1) above we recover the definition of a 
J-diagram [121 Hi], an object which represents a cell in the totally nonnegative part 
of the Grassmannian. 

We will now define a few statistics on permutation tableaux. We define the 
weight wt(T) of a permutation tableau T with k columns to be the total number 
of l's in the filling minus k. (We subtract k since there must be at least k l's in 
a valid filling of a tableau with k columns.) In other words, we are counting how 
many extra l's each column contains beyond the requisite one. 

We define /(T) to be the number of l's in the first row of T. 

We say that an entry in a column of a permutation tableau is restricted if that 
entry is a which lies below some 1. And we say that a row is unrestricted if it does 
not contain a restricted entry. Define u(T) to be the number of unrestricted rows 
of T minus 1. (We subtract 1 since the top row of a tableau is always unrestricted.) 

Figure [3] gives an example of a permutation tableau T with weight 19 — 10 = 9 
and expanse 17, such that u(T) = 3 and /(T) = 5. 

Now we will be interested in the question of enumerating permutation tableaux 
according to their shape, weight, unrestricted rows, and first row. That is, we are 
interested in computing the polynomials F\(q) := q wt ( I 7a~f( I ")0~ u (' r ', where 
the sum ranges over all permutation tableaux T of shape A. Let us also define 
the polynomials F n (q) := ^ r q^i 7 ) ct~f( T ^ 0~ U ( T \ where the sum ranges over all 
permutation tableaux T with expanse n. 

Our main result is the following. 




FIGURE 2. Young diagram and path for A = (2, 1, 0) 
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A = (10,9,9,8,5,2) 



Figure 3. A permutation tableau 

Theorem 3.1. Fix a partition A with expanse n + 1. Lei (t±, . . . , r„) e {0, 1}™ &e 
4>(X). Then 

n 

Fx(q) = ^i(Q( T ^i + (1 - 7-030)^1- 

i=i 

Moreover, the generating function F n+l {q) for all permutation tableaux with ex- 
panse n+1 is W 1 {D 1 + Ei) n Vi. 

An examination of the formula in Theorem 13 . 1 1 reveals that F\(q) is just the sum 
of the entries in the top row of a certain product of n matrices, each one either D\ 
or E\, where this product corresponds via <f> to the shape of the partition A. 

We remark that Theorem [XT] is particularly nice because it is a positive formula. 
In [16j . a recurrence was given for F\(q) (in the case where a = = 1), but not an 
explicit formula. 



Example 3.2. We now illustrate Theorem \3.1\ with an example. Suppose that we 
want to calculate the weight generating function for permutation tableaux of shape 
A = (2, 1). Note that 4>{X) — (0, 1,0). Therefore we need to evaluate the expression 
W\E\D\E{V\. This is equal to a~ 2 (3~ l q + oT 2 + a _1 /3 _1 . And indeed, as shown 
in Figure^ there are three permutation tableaux of shape (2,1), whose statistics 
correspond to the three terms above. Compare our results here to Examvle \2.b\ 
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FIGURE 4. Permutation tableaux of shape (2,1,0) 



We will defer the proof of Theorem 13.11 to Section 2] in order to first explain 
various consequences of the result. 

Theorem 13. II together with Theorem 12.31 implies the following. 

Corollary 3.3. Fix r = (n,...,T„) € {0,1}™, and let A := A(r). (Note that 
expanse(A) = n+1.) The probability of finding the PASEP model in configuration 
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(71, ... , t„) in the steady state is 

Fx(q) 

Here, F\(q) is the weight- generating function for permutation tableaux of shape A. 
Moreover, the partition function Z n for the PASEP model is equal to the generating 
function for all permutation tableaux with expanse n + 1, enumerated according to 
weight, unrestricted rows, and the first row. 

In the case that a = (3 — 1, this corollary together with results from work 
with Steingrimsson i TB] will give us two more combinatorial interpretations of the 
probabilities in the steady state, in terms of excedences and crossings, and descents 
and patterns in permutations. We now recall some definitions and results. 

If 7r £ S n is a permutation, we say that 7r has a weak excedence in position 
i if n{i) > i. Note that there is always a weak excedence in position 1. The 
weak excedence set of tt is the set of all indices i where there is a weak excedence. 
Additionally, we say that 7r has a crossing in positions (i, j) if either j < i < tt(J) < 
7r(i) or < ic(J) < i < j. Crossings were first defined by the first author in [5J. 

Theorem 3.4. |16[ Theorem 7] Let T(k,n,c) be the set of permutation tableaux 
with k rows, n — k columns, and weight c. Let M[k, n, c) be the set of all permu- 
tations 7T £ S n with k weak excedences and c crossings. Then there is a bijection 
$ : T{k,n,c) — > M(k,n, c). Moreover, <1> maps a tableau T whose path (f>(T) has 
south steps in positions S C {1, . . . , n} to a permutation with weak excedence set S . 

Given r = (t%, . . . , r n ) £ {0, 1}™, we define W(t) to be the subset of {1,2,..., n+ 
1} which contains 1 and also contains i + 1 if and only if = 1. Define the excedence 
weight wt(-7r) to be the number of crossings of w, and define F^(q) :— X)tt 9 W * > 
where the sum is over all permutations in S n + 1 with weak excedence set W(t). 

By combining Corollary 13.31 and Theorem 13. 4\ we obtain the following result. 

Corollary 3.5. Fix r = (ti, . . . , r„) G {0, 1}" and suppose that a = j3 = 1. The 
probability of finding the PASEP model in configuration . . . ,r„) in the steady 
state is 

Z n 

Here F^(q) is the generating function which enumerates permutations in S^+i with 
weak excedence set W(t) according to number of crossings. 

This corollary can also be obtained in another way. Namely, in [3], the authors 
provide a lattice path interpretation for the probability of find the PASEP model 
in a particular configuration. If one then applies a certain bijection on lattice paths 
[51 Lemma 9] and then the bijection of Foata and Zeilberger [ID] from lattice paths 
to permutations, one arrives at Corollary 13. 51 in this manner. 

Example 3.6. Suppose that we are interested in r — (0, 1, 0). If we set a = = 1 
and use the results of Example \2.b\ then we find that the probability of finding the 
PASEP model in configuration r in the steady state (for a = (3 = 1) is 2 q i +ik l +iA • 
Note that W(t) — (1,3). Now observe that there are precisely three permutations 
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in Si with weak excedence set (1,3) : (3,1,4,2), and (4, 1, 3, 2), and (2, 1, 4, 3). The 
first of these permutations has one crossing and the second two have none, so the 
generating function for permutations of type r is q + 2, in agreement with Corollary 

E3 

We now introduce some more definitions concerning permutations. If tt G S n is 
a permutation, we say that tt has a descent in position i if tt(i) > 7r(z + 1). The 
descent set of tt is the subset of {1, 2, . . . , n — 1} where 7r has descents. And we say 
that 7r has an occurrence of a generalized pattern of type 2 — 31 if there is a pair 
such that 1 < i < j < n and ir[j + 1) < ir(i) < ir(j). These patterns were 
defined by Babson and Steingrimsson pQ. 

Given r = (n,... , r n ) G {0, l}" -1 , we define D{r) to be the subset of {1, 2, . . . ,n} 
which contains i if and only if r< = 1. We define the descent weight wt'(7r) to be the 
number of crossings of tt, and define F?(q) := <? wt ' (7r) , where the sum is over all 
permutations in SVi+i with descent set D(t). 

We need the following result from [16 , which is also proved in proved in J5[ 
Proposition 6] . 

Theorem 3.7. |16[ Theorem 16]. There is a bijection ^ on permutations in S n , 
which sends a permutation tt with descent set D{t) to a permutation ^{tt) with 
weak excedence set W(t). Moreover, the excedence weight ofty(it) is equal to the 
descent weight of^^ir), i.e. the number of crossings of^{Ti) is equal to the number 
of occurrences of the pattern 2 — 31 in tt. 

By applying this result to Corollary 13. 51 we get the following. 

Corollary 3.8. Fix t = (ti,...,t„) € {0,1}" and suppose a = j3 = 1. The 
probability of finding the PASEP model in configuration (ti,... ,t„) in the steady 
state is 



Here F"(q) is the generating function which enumerates permutations in 5 n +i with 
descent set D(t) according to the number of occurrences of the pattern 2 — 31. 

Now we will show how these results imply a result of the first author [5] . 
First we recall the definition of the polynomials Ek, n (q) which were introduced 
by the second author [TT]. Define 

fc-i 

E k , n (q) = q k - k2 Y,(-mk-i} n q k *- 

It was shown (implicitly) there (and more explicitly in [16) ) that Ek. n (q) enu- 
merates permutation tableaux with k rows and n — k columns according to weight. 
Additionally, it was shown in |17j that at q = —1, 0, 1, Ek tTl (q) specializes to bino- 
mial coefficients, Naryana numbers, and Eulerian numbers. 

In [5], the following connection was made between these polynomials and the 
PASEP model. 
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Theorem 3.9. |5j Let a = (3 = 1. Then in the steady state, the probability that 
the PASEP model with n sites is in a configuration with precisely k particles is: 

Ek+i, n +i{q) 
Z,, 

Since the polynomials Ek, n {q) enumerate permutation tableaux with k rows ac- 
cording to weight, Corollary 13.31 implies Theorem 13.91 

4. Proof of Theorem 13.11 

Recall the following requirement for permutation tableaux. When we read the 
columns of a permutation tableau T from right to left, if in any column we have a 
which lies beneath some 1, then all entries to the left of (which are in the same 
row) must also be O's. 

Also recall that an entry in a column of a permutation tableau is restricted if 
that entry is a which lies below some 1 . And a row is unrestricted if it does not 
contain a restricted entry. 

Proof. We will prove Theorem 13.11 inductively, by finding a precise combinatorial 
interpretation for each of the entries in the top row of a matrix product such as 

n 

Y[{ n Di + (i-Ti)Ei). 

i=l 

Let F{{q) := £ r <7 wt(T) a" /(T) /3~ u(T) = Et <7 wt(T) a" /(T) /3" ,;+1 , where the sum 
ranges over all permutation tableaux T of shape A which have precisely i unre- 
stricted rows. Clearly F\(q) — J2i>i E\(l)- Note that a permutation tableau will 
always have at least one unrestricted row (the top one). 

Now let M\ be the matrix n™=i( r i-Oi + (1 — T~i)E{), where (ti , . . . , r n ) = 0(A). 
We claim that the entry Ma[1, i] in position (1, i) of M\ is F\(q). 

First note that this claim holds when M\ is equal to D\ or E\. M\ is equal 
to Di when A is the partition (0,0). In this case there is a single permutation 
tableaux with two rows and no columns, which has weight 0, and two unrestricted 
rows (hence f(T) = and u(T) = 1). This corresponds to the fact that the top row 
of D\ is (0, f3~ 1 , 0, . . . ). Similarly, M\ is equal to E\ when A is the partition (1). 
In this case there is a single permutation tableaux with one box which is filled with 
a 1 (hence f(T) = 1). This tableaux has weight and one unrestricted row (hence 
u(T) — 0), corresponding to the fact that the top row of E\ is (a" 1 , 0, 0, . . . ). 

Using induction, assume that our claim is true for A with expanse less than 
or equal to n. In other words, we can interpret the ith entry of the top row of 
M\ as a generating function enumerating permutation tableaux of shape A with i 
unrestricted rows, according to weight. Let us now consider how these generating 
functions will change if we instead consider permutation tableaux of shape A', where 
A' is a new shape obtained from A by either adding a new row of length 0, or a new 
column whose length is the number of rows (including rows of length 0) of A. The 
corresponding operation on paths is the following: we take the path corresponding 
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to A and add an additional step from its south-west corner which is either south or 
west. 

It is easy to see how adding a step south to the partition path — i.e. adding an 
empty row to A — will affect the generating functions. Any permutation tableau 
T of the new shape A' will be a permutation tableau of shape A with one addi- 
tional unrestricted row (the last row). Therefore if A' is equal to A union a row 
of length 0, we will have F^(q) = 0, and F^iq) = /3 _1 F|(g) for i > 1. This 
corresponds to the fact that in the matrix product M\Di, the new top row will be 

(O.^MaIMI^-'MaIU],..)- 

It now remains to see how adding a step west to the partition path — i.e. 
adding an extra column of length r to the left-hand-side of A = (Ai, . . . , A r ) 
will affect the generating functions. In this case our new partition A' is equal to 
(Ai + l,...,A r + l). 

Let h a , b (q) denote the polynomial (3 b - a (orV -1 + £°=o q j ( 6 "" +J ')) for 
a < b. We claim that F£,(q) = £ b > a h a . b {q)F b x (q). 

To prove this, it is enough to show the following. Fix a permutation tableau T 
of shape A which has precisely b unrestricted rows. Consider all ways of adding an 
additional (maximal) column C to the left of T, in order to build a new permutation 
tableau T' of shape A' which has precisely a unrestricted rows. Then the generating 
function for these new tableaux according to weight is precisely h a ,b{q)q wt ^ ■ 

To prove this last statement, consider the process of adding the column C to 
T. Since T has b unrestricted rows, there are only b entries in C in which we can 
choose to put either a 1 or 0; all other entries are forced to be 0. Let us number 
those unrestricted entries from top to bottom by c\ to c;,. Since we want our new 
column C to add an additional b — a restricted positions, in our filling of the entries 
ci, . . . , Cb, we must have precisely b — a 0's below a 1. If the top- most 1 in our 
filling of ci, . . . , Cb is in position a, then the i — 1 entries above it are all 0's, and we 
have precisely ( b Zl) ways to choose which entries to make below Cj. The other 
(b — i) — (b — a) = a — i entries must be l's. This particular choice of column 
will therefore contribute the extra weight q a ~ % to the weight of T if i ^ 1 and 
will contribute the weight a~ l q a ~ l - if i = 1. If wc sum over all possible columns 
C which we may add to T, we get the following: arV -1 (^) + ^1=2 Q^ib-i) 
which is equal to oT 1 q a - 1 ( b a Z\) + E"=o 9* ( b ~°j +j ) ■ This completes the proof that 
= Ei> f l a ! b{q)F b (q). And now note that this corresponds to the fact that 
in the matrix product M\E\, the entry in the first row and ath column will be 
Zb> a ha,b(q)Mx[l,b]. 

This now completes our proof that for any partition A, the entry Ma[1,«] in 
position (l,i) of M\ is Ft(q). And since F\(q) — J2i>i the theorem follows. 

□ 

5. Applications of permutation tableaux 

In this section we will give some applications of the connection of permutation 
tableaux to the PASEP. The first application is a partial order on states of the 
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PASEP and some monotonicity results (with respect to that partial order) of prob- 
abilities of observing these states. Note that some of the results below hold for 
general a,/3 (0 < a < 1 and < j3 < 1), and for others we need to assume that 
a = P = l. 

Definition 5.1. Let t,t' E {0,1}" be two states of the PASEP which contain 
exactly k particles. We define the partial order -< by r < t' if and only if A(t) C 
A(r'). 

Figure [5] illustrates this partial order when n = 5 and k = 2, and it shows the 
(unnormalized) probabilities f n (T~) that each of these states occurs, for a = j3 = 1. 



• •ooo 

q 6 +5q 5 +14q 4 +26q 3 +33q +26q+10 



• OiOO 

q 5 +5q 4 +13q 3 +21q 2 +20q+9 



•oo«o 

q 4 +5(f +llq 2 +13q+7 



OftOO 

q 4 +4q 3 +9q 2 +llq+6 



• OOO* 

q 3 +4q 2 +6q+4 



omomo 

q 3 +4q 2 +7q+5 



o#oo# 

q 2 +3q+3 



oo« 

q 2 +3q+3 



>o 



oo«o« 

q+2 



OOOt 

1 



FIGURE 5. A partial order for states of the PASEP 

We now give two simple inequalities relating f n (r) to /„(t') when r -< r'. 

Proposition 5.2. Let a and (5 be general. Suppose that r < r' and let d := 
|A(t')| — |A(t)|. That is, d is the difference between the cardinalities of the Young 
diagrams A(r') and A(r). Then f n { T ') ~ Q d fn(T) is a non-negative polynomial in q, 
a -1 , and (3~ x ■ 

Proof. Observe that any permutation tableau of shape r can be naturally extended 
to a permutation tableau of shape r' by filling in all boxes of A(r') \ A(t) with l's. 
This new permutation tableau has weight q d times the weight of the old one. □ 



14 



SYLVIE CORTEEL AND LAUREN K. WILLIAMS 



Proposition 5.3. Let a = (3 = 1. Suppose that t -< r' . Then f n {r') — fn{T~) is 
a non-negative polynomial in q. In other words, as one moves up the partial order 
-<, the coefficients of f n (r) monotonically increase. 

Propositiori l5.3l is false for general a and j3. For example, /4(1, 1, 0, 0)— /4(1, 0, 1, 0) 
is a Laurent polynomial in q, a, and f3 with some negative coefficients. 

We will defer the proof of Proposition 15.31 to the next section because it is most 
easily proved using Motzkin paths. 

We thank the referee for pointing out that the corollary below follows from our 
previous results. 

Corollary 5.4. Let a — (3 — 1. Suppose that r -< t' . Let d be any integer such 
that < d < | A(t') I — |A(t)|. Then / n (r') — q d fn(i~) is a non-negative polynomial 
in q. 

Proof. This follows from Propositions 15 . 21 and 15.31 Simply choose a saturated chain 
of partitions r = tq -< T\ -< . . . t to = r' such that |A(rfc+i)| — |A(rfc)| = 1 and apply 
either Proposition 15.21 or Proposition 15.31 at each step. □ 

The above results make sense intuitively, since in a model with particles entering 
from the left and leaving to the right, it is more likely that a given particle will be 
further to the left than to the right. Also recall that the probability of hopping left 
(■^zi for some q < 1) is at most the probability of hopping right (^j) 

Another application of permutation tableaux is that these objects allow one to 
read off the main recurrences for the PASEP ([3] Theorem 1]) quite easily. The 
rest of the results in this section hold for general a and /3. 

Fix a partition A = (Ai, . . . , A&), and let r = (n, . . . , r n ) be the unique vector 
in {0, 1}™ such that A = A(r). Choose any corner of A, i.e. the last box of a row 
A; in A such that Aj+i < A^. (Equivalently, a pair (j,j + 1) of entries in r such 
that Tj = 1 and r^+i =0.) As Figure [5] illustrates, there is a simple recurrence for 
Fx(q). 




FIGURE 6. Recurrence for F\(q) 



Explicitly, any valid filling of A is obtained in one of the following ways: 

• inserting a column whose bottom entry is f and whose other entries are 
after the (Ai — l)st column of a valid filling of (Ai — 1, Aa — 1, . . . , \k — 1); 

• adding a f to the end of the ith row of a valid filling of the shape (Ai, A2, • • • , A, 
1, • • Afe); 
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• inserting an all-zero row of length Xi as the ith row in a valid filling of 
(Ai, . . . , Aj-i, Ai+i, . . . , Afc)- 

Note that none of the operations above change the number of 1 's in the first row 
or the number of unrestricted rows. Therefore, for 1 < i < k and Ai > A»+i, wc 
have that 

F\(q) = F { x 1 -i,...,\ i ~i,\, +1 ,...,\ k )+F(\ u ...^ t _ u x i+u ... : \ k )+qF {Xu ... : \ 1 _ u \ l ^i,\ i+u ...^ k )- 

By Corollary 13.31 this recurrence translates into a recurrence for the PASEP 
model. Recall our notation /"( ri ^---' r ") f or the probability of finding the system in 
configuration (ti, . . . , r n ) in the steady state. 

Corollary 5.5. [3J Theoreml] Let a and (3 be general. Then 

/n(Tl, T 2 , . . . ,Tj-i, 1, 0,T 3+2 , . . . ,T„) =/n-l (n , T 2 , . . . , Tj-1,1, Tj+2, ■ ■■ , T„) + 

<lfn(n,T2, ■ ■ .,Tj-l,0,l,Tj+2, ■ ■ ■ ,T„) + 
fn-l(ri, . . . ,Tj-i,0,Tj+2, ■ ■ ■ ,T n ). 

It is also easy to prove the following, by consideration of permutation tableaux. 

Corollary 5.6. [31 Theoreml] Let a and (3 be general. We have that / n (ri, . . . , T n _i, 
jjfnin, ■ ■ -,T n -i) and also /„(0,ti,t 2 , . • . ,r n ) = ~/„(ti,t 2 , . • . ,r n ). 

Proof. This also follows from Corollary 13.31 The first equality holds because the 
generating function for permutation tableaux of shape (Ai, A2, . . . , A„_i, 0) is clearly 
equal to the generating function for permutation tableaux of shape (Ai, A2, . . . , \ n -i) 
except that there is one additional unrestricted row. The second equality holds be- 
cause the generating function for permutation tableaux of shape (Ai + 1, A2, . . . , A„) 
is equal to the generating function for permutation tableaux of shape (Ai, A2, . . . , A„) 
except that there is one additional entry in the first row which is a 1. (Recall that 
each column of a permutation tableaux is required to contain exactly one 1 and 
that the weight of a column is equal to the number of l's it contains beyond that 
requisite one.) □ 



6. Connection with Motzkin paths 

We now give another solution to the matrix ansatz, which is essentially the one 
described in Derrida et al [7] , and which has an interpretation in terms of Motzkin 
paths. 

Note that in this section we restrict to the case a — (3 = 1 . These results can be 
extended to the more general case, but we restrict ourselves to this case since our 
main purpose here is to obtain the result needed for the monotonicity result of the 
previous section. 

Let Dq be the (infinite) upper triangular matrix (dij) such that di^ = [i] where 
[i] is the g-analog of i, = [i + 1], and dij = for j : ^ i, i + 1. 
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[1] [2] ... \ 

[2] [3] ... 

[3] [4] ... . 

[4] ... 

V i ) 

Let Eq be the (infinite) lower triangular matrix (e^j) such that = [i], e^+i^ = 
[i], and e^j = for j ^ i, i — 1. 
That is, Eq is the matrix 



1 











1 


[2] 











[2] 


[3] 











[3] 


W 



V ; ; i ; j 

Let Wo be the (row) vector (1, 0, 0, ... ) and Vo be the (column) vector (1, 0,0,...). 
It is now easy to check that D Q E — qE Do = D a + E , DV = Vo, and W Q E = W . 

We will now give an interpretation of the steady states of the PASEP in terms 
of bicolored Motzkin paths. This result is very similar to a result obtained in Brak 
et al [3], via a combinatorial derivation. Note also that in [2], Brak and Essam 
considered the case q = of the PASEP and gave multiple interpretations for the 
steady states in terms of various weighted lattice paths. 

We define a bicolored Motzkin path of length n to be a sequence of steps in the 
plane c = (c±, . . . , c n ) such that Cj £ {N, S, E, E} for 1 < i < n, which starts and 
ends at height 0, and always stays at or above height 0. That is, if we define the 
height at the zth step to be hi = {j < i Cj = N} — {j < i Cj — S}, then hi > 
for 1 < i < n, and furthermore, h n+ \ = 0. 

We now assign a weight wt(p) to a bicolored Motzkin path p as follows: step Ci 
is assigned weight [h + 1] if Ci ends at height h, and the weight of c is defined to be 
the product of the weights of all the c, 's. 

Fix a sequence r = (ti, . . . , r n ) £ {0, 1}". We say that a bicolored Motzkin path 
c has type t if a — N or E whenever r, = 1 and c, — S or E whenever Tj = 0. 

We now define F"' (q) to be the generating function for all bicolored Motzkin 
paths of type r: that is, F"'(q) :— X) P '? wt( ' P ' > where the sum is over all bicolored 
Motzkin paths of type r. 

The main result of this section is the following. 

Proposition 6.1. Choose a sequence r = (n, . . . ,r„) £ {0, 1}". Then 

n 

F';\ q ) = WoiHinDo + (1 - n)E ))Vo 

i=l 

Let us make the convention in this section that rows and columns of our matrices 
are indexed by non-negative integers (including 0). Then Proposition 16.11 says that 



That is, Do is the matrix 
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the generating function F"'(q) is equal to the entry in the Oth row and 0th column 
of the matrix product FKLi O 7 *-^ + (1 — t^Eq). 

Proof. Using the definition of matrix multiplication, we expand n™=i( T i-^'o + (1 — 
ti)Eq) as the sum of terms of the form At[l,ii]i4.2[ii, £2] ■ ■ -A n [i n -i, 1] where the 
Aj's are matrices Dq or Eq, depending on the sequence r. It is obvious that these 
terms correspond to the Motzkin paths of type r. □ 

As before, we can use Proposition 16.11 together with Theorem 12.31 to deduce the 
following. 

Corollary 6.2. Fix t — [t\, . . . ,T n ) E {0,1}". The probability of finding the 
PASEP model in configuration (ti, . . . , r n ) in the steady state is 

Here F"'(q) is the generating function for bicolored Motzkin paths of type r. 
We now use this result to prove Proposition 15.31 from the previous section: 

Proof. It is sufficient to consider two configurations r and r' whose corresponding 
diagrams A(r) and A(r') differ by one box. That is, there exists some £ such that: 

• Ti = 0, T l+ l = 1 

•t! = 1,t! +1 = 

• Tj = Tj for j =/= i, i + 1. 

We now observe that any bicolored Motzkin path p of type r (with weight w) can 
be mapped to a bicolored Motzkin path p' of type t' whose weight w' is coefficient- 
wise greater than w. 

If p = (ci, . . . , c„) then we let p' := (ci, . . . , Cj_i, c i+1 , Q, c l+2 , . . . , c„). That is, 
p' is the path obtained from p by switching the zth and i + 1st steps. It is clear 
that whether we switch SN steps for NS, or SE for ES, or EN for NE, or EE 
for EE, the resulting path is a valid Motzkin path. The weight w' will differ from 
w by (respectively): replacing a factor [h — l][h] with [h + l][h]; replacing a factor 
[h — l] 2 with [h][h — 1]; replacing a factor [h][h + 1] with [h + l] 2 ; nothing. Using 
the fact that [n + 1] — [n] is nonnegative, all of the differences [h + 1] [h] — [h—1] [h] , 
[h] [h—1] — [h—1] 2 , [h+ 1] 2 — + are nonnegative; therefore in all cases, w' — w 
will be a nonnegative polynomial. By Corollary 16. 21 we are done. □ 

Note that one can use the same argument to give another proof of Proposition 
15.21 in the case that a = j3 = 1. We repeat the same argument but now w' differs 
from w by: replacing q[h — l][h] with [h + l][h]; replacing q[h — l] 2 with [h][h — 1]; 
and replacing q[h] [h + 1] with [h + l] 2 . Then all differences of such polynomials will 
be nonnegative, since [n + 1] — q[n] is nonnegative. 

Remark 6.3. We now summarize our combinatorial interpretations for the steady 
state probabilities of the PASEP model; these are given by our solutions (Dj,Ej) 
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to the matrix ansatz (where a = (3 = 1), for j = 0,1. Note that since 

n 

W 3 ([[(nD 3 + (l - n )E 3 ))V 3 

i=l 

describes the steady state probability of the PASEP model for both j = 0, 1, these 
two formulas must be equal. In particular, both of these are formulas for all of the 
following: 

• F\(r)(q), the weight-generating function for permutation tableaux of shape 
A(r). 

• F^iq), the generating function for permutations of excedence type t, enu- 
merated according to crossings. 

• F"(q), the generating function for permutations of descent typer, enumer- 
ated according to occurrences of 2 — 31. 

• F"'(q), the generating function for weighted bicolored Motzkin paths. 

Moreover, the product Wj(Dj + Ej) n Vj is the weight-generating function for all 
of the following: 

• permutation tableaux with expanse n + 1 . 

• permutations in S n +i, enumerated according to crossings. 

• permutations in S n+ i, enumerated according to occurrences of 2 — 31. 

• weighted bicolored Motzkin paths of length n. 
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